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Abstract 

The temperature and pressure effects on the resistivity of the manganese ox- 
ides are studied analytically via the Kondo lattice model with a strong Hund's 
coupling. We obtain analytically the single-particle density of states on the 
Bethe lattice in the large connectivity limit using the analytical variant of the 
dynamic Lanczos method. From the density of states for the doped system, 
we obtain resistivity and show resistivity drop when temperature crosses the 
magnetic transition point. We also demonstrate the effect of pressure both 

below and above the transition temperature. 
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I. INTRODUCTION 



Though manganese oxides had been studied extensively during 50 's and 60 's, [[[] the 
hole-doped manganese oxides La\- x A x MnOz(A = Ca, Ba, Pb, or Sr) have recently driven a 
considerable attention due to not only the phenomenon of giant magnetoresistance (GMR) 
but also various interesting physical properties such as doping induced metal-insulator 
transition, antiferromagnetic-ferromagnetic transition, and the structural phase transition 

These phenomena have been partly explained by the double exchange mechanism. 
The ferromagnetic Kondo lattice model (KLM) due to strong Hund's coupling |]] has been 
proposed to describe the properties of manganese oxides. Recently, a mean field theory, 
which is valid in infinite dimensions, Q has been applied to the KLM in the limit of infinite 
magnitude of local spin. || However, the large resistance drop measured in experiment 
has not been shown successfully. Millis et al. have recently argued that the strong electron- 
phonon coupling is necessary to explain the colossal magnetoresistance ||. a new experiment 
[10|1 shows that there is no primary lattice distortion anomaly at the magnetic transition 



temperature. 

Leaving aside the question how much the polaron physics affects the properties of man- 
ganese oxides, we are interested in how well the KLM with a strong Hund's coupling de- 
scribes the resistivity drop upon variation of temperature [II] and pressure |)T2]-fH| in the 
manganese oxides. We feel that proper calculation of the resistivity in this model is necessary 
to understand how much the electron-phonon interaction is important for the explanation 
of the experiments. Even though magnetovolume contraction accompanies across 

the magnetic transition, we do not consider this effect explicitly in this work. We suppose 
that the effect of pressure is to enhance not only hopping between conduction electrons but 
also scattering rate. These two effects appear differently in ferromagnetic and paramagnetic 
states. The purpose of this work is to show quantitatively how much the resistivity drops 
due to magnetic ordering and pressure via temperature and transfer integral, respectively. 

Our work starts from the single-particle density of states (DOS) which has not been 



2 



previously obtained through a rigorous calculation for the KLM. We first obtain analytically 
the DOS of the KLM on the Bethe lattice in the large connectivity. The method we use 
here, the so called analytic version of the dynamic Lanczos method, has proved useful in 



describing the optical transport properties of V2O3 [17|] and the doped cuprates |]l8| . In a 
previous work ||T^|, analytical calculation of DOS through this method has given transparent 
understanding of the metal-insulator transition in the large-dimensional half-filled Hubbard 
model . 

This paper is composed as follows. In Sec. II, we briefly introduce the formalism, and 
show the single-particle DOS in Sec. III. We present the resistivity depending on temperature 
and pressure in Sec. IV, and finally give conclusion in Sec. V. 

II. FORMALISM 



The single-particle DOS, N a (u), is given by |I9| , |22 | 

N a (u) = — lim £ Re / ({c ja (t), c) a })e^dt 

7TiV T7-+0+ j JO 

= - lim Re a (z)| z =_^+»? ; (!) 

7T 77— *0+ 

where the angular brackets denote the thermal average and a (z), which is the Laplace 
transform of ({cj a (t), Cj a }). One notes that a (z) corresponds to the on-site Green function 
and is the key quantity giving the DOS. 

One can observe that ao(z) is just the projection of Cj a (z) onto Cj a in a Liouville space 
defined by the inner product (A, B) = ({A,B^}), where A and B are operators of the 
Liouville space, and B^ is the adjoint of B. This projection can be obtained most easily in 
the orthogonalized space. Therefore, we first construct an appropriate orthogonal space. We 
choose the first vector / as fo = Cj a and the second using the Gram-Schmidt process with 
a linearly independent vector iLf where L is the Liouville operator. The rest orthogonal 
vectors are obtained via a similar process which may be expressed by the recurrence relation, 

= iLf v - a v f v + A v f v -i, (2) 



where /_i = and A u = (f v ,fv)/{fv-i,fv-i) for v > 1. Use of the relation (iLA,B) = 
—(A,iLB) has been made. Applying Eq. (2) to f (t) = J2^Lo a u(t)fu, the projection a (z) 
is represented by the infinite continued fraction [ET]-|2H] 



= ; - i 1 a t: — ' (3) 

Z - "0 + z _ ai+ A 2 



where a v = (iX/„, /„)/ (/,/, /„)■ 

For a nontrivial system, however, the fractional forms of a v and A u require expansions in 
performing the Gram-Schmidt process analytically, and allow us an orthogonal set satisfying 
Eq. (2) exactly only in an asymptotic regime, e.g., strongly or weakly correlated regime. 
Since a v provides with energy standard and can be controlled by chemical potential, a v 
usually has a simple form. What makes things hard is A u which is the ratio of norms of 
consecutive bases. Since norm is usually given by a polynomial composed of positive terms of 
an expansion parameter, the ratio of leading terms gives quite good approximation when the 
terms of each norm show a similar behavior. This is actually the case we meet in practice. 
This maybe happens since we get the bases through a simple recurrence relation. We argue 
that the valid region of the asymptotic dynamics extends to the regime where the expansion 



is allowed. One can recognize it from the result. One of us and Lee p6|j27[] have shown the 
dynamics of electron gas in the asymptotic regimes such as small momentum transfer and 
large momentum transfer limit. They have shown that the former is good for the regime 



k < kp [p6 |, where k and kp are momentum transfer and Fermi momentum, respectively 



and the latter agrees well with experiment for the momentum transfer k ~ 2kp [[27 . 



III. SINGLE-PARTICLE DENSITY OF STATE 

We now apply this formalism to the KLM which is written as 

H = - 1 >j { - r ir r j.n + h.c.) - JJ2 s i " °h (4) 
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where angular bracket means nearest neighbor. The first term represents hopping processes 
of the itinerant electrons in e g orbital and the second represents ferromagnetic Kondo cou- 
pling between localized spin (S = 3/2) and itinerant electron spin on Mn 3+ ion coming from 
t 2g and e g orbitals, respectively. The electron spin Oi can be represented in terms of fermion 
operators. Therefore we have the following Hamiltonian. 

H = - S *«( c Jt c j> + h - c -) - \ Yl s iA\ c a + s tAfi\ + s i ( c ! T c *t - c \i c a)\ ( 5 ) 

where S = (S x , S v , S z ) is a classical spin. 

The key part of this work is to construct the orthogonal bases f^s along with a u and 
A v . Using the recurrence relation (2) with f = Cj<\ and considering J » Uj for expansion, 
we obtain the orthogonal bases by taking the leading terms as follows: 

1v 2v 2v 

fiu = J 5]j ' ' ' ^> q tkj ' ' ' tpg Tfc • • • T p C q ^ 

2^ 2v 2v 



fiv+l = i ~ £j • • • S ? t k j ■ ■ ■ t pq T k ■ ■ ■ T p (ASgCq1 + S q c ql ) 



2v+l 2v+l 2v 



fav+2 — —J 2U+1 £j • • • S' r tkj • ■ - t qr Tk - • ■ T p [(Sg + S*)c r <[ + (S q + S r )c r J 



fiv+z — — * — 2 — ^jf ' ' ' ^kj ' ' 'tqr Tfc ' ' ' Tp 

x[(S z q + S z )AS z c r ^ + (Sq- + ,S r 7)5 r + c rT + (S z q S z r - S q S z )c ri 

-(s- + S;)(S z )c rl ] (6) 

where v > 0, AS 2 = S 2 - (S 2 ), = (S^ 2 ) = (S 2 ), and T k = S k + S,. We consider only 
nearest neighbor hopping so that j is the nearest neighbor of k, I is the nearest neighbor 
of j, and the prime in summation means k 7^ j 7^ • • • 7^ r. The inner operation of {T^} 
product is defined as the scalar product only with itself. This operation is valid for classical 
spins. We obtain the inner products as follows 

J 2 



(hJi) = (" T [(&s z k ) 2 + s k -sz\) 
j 2 



= ^((sl) - (s k ) 2 ) 



(/ 2 , h) = (J%tl 3 \{si + s;f + (s k + s-)(si + s})\) 

= ^s;<-((s fe + s J ) 2 ) 

= J%t%(Tl) 

(/3,/3) = ^s;4(t2)((s?)-(s,) 2 ) 
(/4,/4) = j 4 s;s^,4(t^ 2 ) 



(7) 



These bases (6) satisfy the orthogonal condition as follow. 

2v 1v 1v 

= {^^~^A^ q n^ P \( s i - (s z )){c qh ci,} + s;{ Cql ,c^}}) 

2v 2v 2u 

^—^ ^ 



^■■■■Kti J ---t' pq (Ti)...(Tp((S z q -(S z })} = (8) 

and so on. We take the approximation neglecting the spatial correlations between the local- 
ized spins of Mn 3+ , that is, (S fc -Sj) = (S fc ) • (Sj) = (S) 2 for k ^ j. where (S) is independent 
of sites. Therefore (T 2 ) = 2(S 2 ) +2(S) 2 . Note that this mean field approximation (MFA) is 
valid in the limit of large connectivity. Within this approximation, one obtains a v and A u 
as follows: 



&2u 




ia, 




- -4<s> 


= — ia, 


^2v+l 




- (S) 2 ) 


A.2v+2 


_ 4t MS 2 ) 

*((S 2 > 


+ (S) 2 ) 
- (S) 2 ) 



= A G , 

= Ae, (9) 

for v > 0, where t* = and q is the coordination number which is taken to be very large. 

Now the DOS of the up-spin conduction electron at site j is given by Eqs. (1), (3), and 
(9), i.e., 



i V 4 ^ 2 - « 2 ) - ( A ° - A * - - « 2 )) 2 

Nt(uj) = — — — (10) 
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This DOS shows two separated bands of width W = J a 2 + (yAo + v^e) 2 _ 



a 2 + (y/K^ — \TK~ e ) 2 , which is 4£* when (S) = 0, i.e., in paramagnetic state. Eqs. (9) 
and (10) show that the DOS depends on the magnetic properties of the system through (S) 
and (S 2 ) of the localized spin. 

Below the transition temperature T c , the system is ferromagnetic and the expectation 
value of (S) can be obtained from the expermental data and (S 2 ) is nothing but S(S + 1). 
Above T c , where the system is paramagnetic, (S) becomes while (S 2 ) remains S(S + 1). 

Fig. 1 shows the DOS for the Hund's coupling strength J/W = 4 and a doping concen- 
tration x = 0.3 using (S) = 0, 0.4, and 0.8 with dotted, dashed, and solid lines, respectively. 
The vertical lines indicate positions of the chemical potential. We only show the lower band 
since the upper band is nothing but point reflection of lower band about zero frequency, 
and does not affect resistivity within the range of temperature in which we are interested. 
The DOS clearly shows the shift of the spectral weight to lower energy side as local spins 
are ordered. Magnetic ordering enhances up-spin spectral weight and diminishes down-spin 
spectral weight remarkably. Note that the double exchange model leads us to the rigid band 
picture and the position of chemical potential is determined only by the amount of filling 
which results in the change of resistivity upon doping concentration, though our DOS, Eq. 
(10), does not contain the carrier concentration (n). [|I{J Our DOS satisfies sum rule exactly. 

IV. RESISTIVITY 

The physical quantity we want to examine in this work is resistivity affected by ordering 
of magnetic moments and pressure. We show specifically in what follows how magnetic 
ordering and pressure affect resistivity via obtaining optical conductivity. 

Optical conductivity in infinite dimensions is obtained using the following formula [30]. 

a(u) = a 1*J dw' J deA (e)A(e, u')A(e, J + u) f ^'' ~ f J W + ^ = a I(u), (11) 



where A (e) = ^2 — p-/t*7r for the Bethe lattice, f{u) is the Fermi distribution function, 



and do = ne l£y N where a, N, V are lattice constant, number of lattice sites, and volume, 
respectively. 

One can get A(e,u) by calculating the self-energy £(a>) using the relation T^ a (uj) = 



to 



2 G a {w) 



applicable to the Bethe lattice |T7| , |3T| and G(u) = —ia (—iuj). The dc 



resistivity is obtained by taking the zero-frequency limit for the inverse optical conductivity, 



i.e., 



p{T) = lim = 1 (12) 

The unit of conductivity is provided by <7o which, in three dimensions, approximately gives 
Mott's minimum metallic conductivity auott ~ 10 3 (ficm) -1 . |8|,|30| However the absolute 



value of the resistivity should depend on the microscopic parameters such as lattice constant 
in cr , so one cannot determine a precisely. It usually varies of 0(1) depending on sample. 

Fig. 2 shows temperature dependences of the resistivity for J/W = 4 with x = 0.3 
and magnetization around transition temperature. The magnetization is obtained as M = 
g[(S) + |((^t) ~ ( n l))]i where g is the g-factor. The temperature dependence of (S) has been 



conjectured as (S) = y (S 2 )(l — |r) from the mean-field type result, i.e., (S) oc yjl — |r. 
This magnetic ordering curve shown in Fig. 2 is remarkablely similar to the experimental 
data of magnetization [0. It shows that the larger magnetic ordering induces the larger 



resistivity drop. Since it has been known that La\^ x Sr x MnO^ has a bandwidth W ~ leV 
HH , we set t* = 0.25eV so that the bandwidth W = 4t* = leV in the paramagnetic state, 



and J = 4W = 4eV which is known as a reasonable value of Hund's coupling. || The 
transition temperature is set T c ° = 0.13t* = 377 K to compare with the experiment JlO . 



We try to understand resistivity drop more physically through the changes of scattering 
rate and bandwidth versus local spin fluctuation. In Fig. 3, we show the change of scattering 
rate, 1/r-f, and bandwidth across T c to see what governs the resistivity drop. The scattering 
rate is given by l/r-j- = 2ImTi^{u = fi), where \x is chemical potential. In Eq. (9) and (10), 
we see that the DOS has a local spin fluctuation term ((AS) 2 ). The local spin fluctuation 
decreases as magnetic ordering increases, which results in broadening of bandwidth and 
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decreasing of scattering rate of up-spin conduction electron which is parallel to the localized 
spin. The scattering rate of down-spin conduction electron increases quite rapidly as local 
spins are ordered. Fig. 3 demonstrates that more spin-disorder scattering occurs in the 
paramagnetic state, and the spin-disorder strongly affects the scattering rate governing the 
resistivity. 

We show pressure effect by changing the transfer integral £* for a fixed Hund's coupling 
strength J. Though the pressure affects both J and bandwidth W, its effect is much larger 
for W than for J. Therefore we make J fix and W increase as pressure increases. The 
proportionality T c cx W is valid under MFA |6|f|, so that T c is taken at A1AK for W = l.leV 
compared with 377 K for W = l.OeV. Fig. 4 clearly shows that no pressure effect appears 
in the resistivity above T c . The origin of this result can be understood as follows. We obtain 
the effective mass using m* = 1 — ^ s \ U=/M . The effective mass decreases by increasing 
t*, while the scattering rate increases nearly the same amount by increasing t*, so the 
effect of change of t* in the resistivity disappears in the paramagnetic state, if we consider 
p oc m*/r with assuming constant number of charge carriers. Therefore we have to consider 
the lattice constant appearing in cr , which will be changed under pressure, to explain the 
experimental results showing remarkable resistivity drop even above T c under pressure. The 
lattice dynamics may play an important role in the paramagnetic regime near T c . || 



V. CONCLUSION 

We study the temperature and pressure effects on the resistivity of the hole doped man- 
ganese oxides by calculating the DOS and optical conductivity for the KLM with strong 
Hund's coupling using the Lanczos continued fraction formalism along with MFA. These 
two effects appear in both DOS and resistivity, and the results show that the smaller spin 
disorder induces the smaller resistivity. Our result also shows that there is no resistivity 
drop in paramagnetic state, above T c , under pressure, which has a clear discrepancy with 
many experimental data |T2]-irjJ]. It demonstrates that either double exchange model losts 
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the main physics in the GMR phenomenon or treating the double exchange model within 
MFA yields too simple physics, which shows that only local spin fluctuation can change resis- 
tivity. Above and near T c , however, the magnetovolume effect may cause further increasing 
of resistivity due to decreasing of lattice constant in o"o of Eq. (12). This effect should be 
included to improve the resistivity drop just above T c and to show the pressure effect shown 
in experiments []T2 -jl^Jl , to some extent. Thus present results indicate that lattice distortion 



may be crucial as far as the resistivity in paramagnet is concerned and polaron may play an 
important role in this sense. 
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FIGURES 

FIG. 1. DOS when J/W = 4 for paramagnetic state (a) with dotted line and ferromagnetic 
state (b) with dashed line and (c) with solid line. We use (S) = 0.4 (b) and 0.8 (c). The vertical 
lines indicate the position of the chemical potential for x = 0.3. 

FIG. 2. Resistivity and magnetization for J/W = 4 and x = 0.3 are shown around the transition 
temperature. Magnetization has been obtained as g[(S) + \{{n-\) — (nj.))], where g is the ^-factor. 

FIG. 3. Resistivity (a), bandwidth (b) and scattering rate (c) in units of 10 — 3 S1 cm, 2t*, and 
g sec." 1 respectively v.s. local spin fluctuation are drawn. 

FIG. 4. Resistivities for W = 1.0eF(solid line) and l.leV(dashed line) with J = 4eV, x = 0.3 
are shown around the transition temperature, where T° denotes the transition temperature for 
W = l.OeV 
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